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Abstract. In this paper we study homomorphisms of Probabilistic Regula- 
tory Gene Networks (PRN) introduced in [2- The model PRN is a natural 
generalization of the Probabilistic Boolean Networks (PBN), introduced by 
I. Shmulevich, E. Dougherty, and W. Zhang in 1141 . that has been using to 
describe genetic networks and has therapeutic applications, see 1151 . In this 
paper, our main objectives are to apply the concept of homomorphism and 
e-homomorphism of probabilistic regulatory networks to the dynamic of the 
networks. The meaning of e is that these homomorphic networks have similar 
distributions and the distance between the distributions is upper bounded by 
€. Additionally, we prove that the class of PRN together with the homomor- 
phisms form a category with products and coproducts. Projections are special 
homomorphisms, and they always induce invariant subnetworks that contain 
all the cycles and steady states in the network. Here, it is proved that the 
e-homomorphism for < e < 1 produce simultaneous Markov Chains in both 
networks, that permit to introduce the concept of e-isomorphism of Markov 
Chains, and similar networks. 



Introduction 

Genes can be understanding in their complexity behavior using models according 
with their discrete or continuous action. Developing computational tools permits 
describe gene functions and understand the mechanism of regulation |H1 E] ■ This 
understanding will have a significant impact on the development of techniques for 
drugs testing and therapeutic intervention for treating human rfzseases[7l 1131 [T5] . 

We focus our attention in the discrete structure of genetic regulatory networks, 
instead of, its dual moving continuo-discrete. Probabilistic Gene Regulatory Net- 
work(PRgN) is a natural generalizations of the model Probabilistic Boolean Net- 
work (PBN), introduced by I. Shmulevich, E. Dougherty, and W. Zhang in jl4| . 
The mathematical background of the model PgRN, is introduced here, for simplic- 
ity we work with functions defined over a set X to itself, with probabilities assigned 
to these functions. X is a set of states of genes, for example X = {0, 1}", if our 
network is a Boolean network. Working in this way, we can observe the dynamic 
of the network indeed focus our attention in the description of functions. The set 
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X can be a subset of {0, 1}", and we can extend some classical ideas to regulatory 
network, such as invariant subnetworks, automorphisms group, etc. In particular 
if X is a vector space over a finite field, the functions are lineal functions, then 
we can use linear algebra to describe the state space. Mapping are important in 
the study of networks, because they permit to recognize subnetworks, in particu- 
lar determine when two networks are similar or equivalent. Special mappings are 
homomorphisms and e-homomorphisms, we use both to describe subnetworks and 
similar networks. An homomorphism transform a network to another in such a 
way the discrete structure giving by the first network can lives in part of the other 
one, or these two networks are very similar but no equals, in particular in the 
probabilistic way. An e-homomorphism is the same but with the condition that 
the probability distributions of the networks are close, and we use a preestablishes 
< e < 1 as a distance between the probabilities. For concept of homomorphism 
of discrete dynamical systems see |51 HI 112). 

1. Preliminaries concepts: Finite dynamical systems, probabilistic 
Boolean networks and Probabilistic Regulatory Networks 

Two finite dynamical systems (AT, /) and (Y, g) are isomorphicas(or equivalents) 
if there exists a bijection (f) : X such that (j> o f — g o cj)^ (or/ = ° .9 o <!>)■ 
If (j) is not a bijection map then cf) is an homomorphism. 

If y C a: is such that f{Y) C Y then (Y, f\y) is a sub-FDS of (AT, /), where /|y 
is the map restricted to Y. There exists naturally an injective morphism from Y to 
X called inclusion and denoted by l. The state space of a FDS (AT, /), is a digraph 
whit vertices the set X, and with an arrow from u to t; if f{u) = v. 

For example, the FDSs X = ({0,l}2,/i(x, 2/) = (rry, y)), and F = ({0, Ip, ^(x, y) 
(x, {x + l)y)) are isomorphics, because their state spaces are isomorphics. 

(1.0) -./i (0,0) (0,1) 

(1,1) 

(1.1) -^f- (1,0) (0,0) 

(0,1) 

In fact, the isomorphism : {0,1}^ {0,1}^ is the bijection 0(1,0) = (1,1), 
0(0,0) = (1,0), 0(0, 1) = (0,0), and 0(1, 1) = (0, 1). The following is an example 
of homomorphism (inclusion) with Z = {{(0, 0), (1, 0)}, /i}. 

(1,0) (1,0) (0,1) 

i/i ifi r^f' 

(0,0) O (0,0) O (1,1) 

A Probabilistic Boolean Network A — (V, F, C) is defined by the following sort 
(type) of objects ?]: a set of nodes (genes) V — {xi, . . . , a;„}, Xi G {0,1}, 
for aU i] a family F = {Fi, i^2, • ■ ■ , F„} of ordered sets Fi = {fi\f2\ 



of Boolean functions /j*'' : {0, 1}" — > {0, 1}, for all j called predictors; and a list 
C = (Ci, . . . , C„), Ci = {c^i \ . . . , c^(^j)}, of selection probabilities. The selection 
probability that the function /j'' is used for the vertex i is Cj*' — Pr{f^'^^ = fj^^}- 
The dynamic of the PBN is given by a vector of functions ffe = ifll\fj^\ • ■ • , /fc"^) 
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for 1 < fci < and /^,*'' G Fi, where k — [fci, . . . , /c„], 1 < ki < The map fk ■ 
{0, 1}" {0, 1}" acts as a transition function. Each variable Xi G {0, 1}" represents 
the state of the vertex i. All functions are updated synchronously. At every time 
step, one of the functions is selected randomly from the set Fi according to a 
predefined probability distribution. The selection probability that the transition 
function = ifil\fj.^J, ■ • ■ , /l"^) is used to go from the state u G {0, 1} to another 
state ffc(ii) = V £ {0, 1}" is given by 

n 
i=l 

The dynamical transition structure of a PBN can be described by a Markov chain 
with fixed transition probabilities. There are two digraphs structures associated 
with a PBN: the low-level digraph F, consisting of genes functions essentiality 
relations; and the high-level digraph which consists of the states of the system and 
the transitions between states. The matrix T associated to the high level digraph 
formed by placing p{u,v) in row u and column u, where m, u G {0, 1}" is called 
the transition probability matrix or chain matrix, p{u, v) = X]fk|fk(u)=f ''ffc' 

1.1. Probabilistic Regulatory Gene Networks. A Probabilistic Gene Regu- 
latory Network (PRN) is a triple X = (X, F, C) where X is a finite set and 
F = {/i, . . . , /„} is a set of functions from X into itself, with a list C = (ci, . . . , c„) 
of selection probabilities, where q = p{fi), HI Q We associate with each PRN a 
weighted digraph, whose vertices are the elements of X, and ii u,v G X, there 
is an arrow going from m to w for each function fi such that fi{u) = u, and the 
probability Ci is assigned to this arrow. This weighted digraph will be called the 
state space of X. In this paper, we use the notation PRN for one or more networks. 



Example 1.1. If X ^ {0,1}\ F = {fi{x,y) = (x, y), ^(x, y) = {x,0), 

f3{x,y) = {l,y)j4ix,y) = (1,0)}; and C = {.46, .21, .22, .11}, the state space of 

X = (AT, F^ C) is the following: 



•67 , 

O (0,0) 
0(1,0)^(1,1)0-68 



(0,1) o-« 
/.111-'' T = 



67 





.33 





21 


.46 


.11 


.22 








1 











.32 


.68 



1.2. Homomorphisms and e-homomorphisms of PRN. If C is a set of se- 
lection probabilities we denote by x the characteristic function over C. That 
is X : C" U {0} ^ {0, 1} such that x(c) = 1, if c 7^ and x(0) = 0. Let 
Xi^iXi,F= (/OLi, C) and ^"2 = (A2, G = (5j),"Li, ^) be two PRN. 



Definition 1.2 ( Homomorphisms of PRN). A map : Ai — > A2 is an homo- 
morphism from Xi to X2, if for all fi there exists a gj, such that for all u, v in 

(1) <P° fi = 9j o 0; and (2) x{dg.{(l){u),(l){v))) > x{cfju,v)). 

Ai ^ Ai 

0i 10 
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(3) Condition for e-Homom,orphism: The distributions of probabilities following the 
homomorphism are enough close. An e- homomorphism is an homomorphism that 
satisfies the condition, for all i, j, max\p{ui, Uj) —p{4>{ui), (j){uj))\ < e, where e > 
is a real number that we previously determine for the applications. 

If <j) : Xx ^ X2 is a bijective map, and for all ft, gj, u, and v in X\; dg^ {(f){u), (l>{v)) = 
Cf^{u,v), then (j) is an isomorphism. 

Example 1.3. 

If A" = (X; F; C) is the PRN in Example 1, and Xi = (X; F' = {/i, /2, /s}; C" = 
{.47, .28, .25}) is a new PRN over the same set X with different probabihties and 
only three functions. 

Xi X 

.47 46 

o (0, 0) ^-28 (0, 1) ^ (0, 0) ^-21 (0, 1) o 

.25 I I .25 .33 I / I .22 

Qi (1, 0) ^ (1, 1) 0-7' (1, 0) 4i (1, 1) C-*^* 



Ti = 



.75 





.25 







.67 





.33 





.28 


.47 





.25 


,T = 


.21 


.46 


.11 


.22 








1 











1 











.28 


.72 










.32 


.68 



The homomorphism (A : A*! — > A" is a bijective map, (j>{x) = x, over the set of states, 
but an inclusion over the set of arrows, because the arrow going from (0, 1) to (1, 1) 
in X doesn't appear in Xi. The first condition for homomorphism is obvious. The 
condition (2) holds, because the inclusion of arrows. The two transition matrices 
are connected by this inclusion, since if the place ij in the first matrix ^ then this 
place is 7^ in the second network too. The two PRN are not isomorphics because 
the probabilities are not equals. In order to determine e for the homomorphism, 
we use the transition matrices. In this example e = .11. 



Ti - r = 



08 





-.08 





07 


.01 


-.11 


.03 




















-.04 


.04 



If the homomorphism is a bijective map like here, the transition matrices Ti and 
T2 have the same order, and X^"=i(T'i — T2)ij = 0, for j = l,n 



2. Applications to Markov Chains, c-Similar Networks 

Two PRN are e-similar if there exists a bijective homomorphism </> between them, 
such that is also an homomorphism. Observe that </> and have the same e. 

Example 2.1. 

The following networks arc .005-homomorphics with the injcctivc homeomor- 
phism : {0,1}^ {0,1}^, giving by (t>{x,y) = {x,y,l). We can see in this 
example, that the network Xi is an e-subnetwork of X2, the special subnetwork 
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4>{Xi) is an invariant subnetwork of A2. The networks, and the transition matrices 
are the following: 



00 



.111, .451 



.549 



10 

.445 

i 

01 

.013.662 

Ti 
11 



000 



.549 

A51 



Xo 



100 

.995 
i 

010 

.378 
O 



.444 

o 



.338 

o 



.987 

o 



.005 



.622 



Ti 



001 
110 

Q.998 






.549 


.451 








.338 





.662 


111 


.445 


AM 








.013 





.987 


.456, 


.113 


.448 
101 





.544 



.002 



.439 
i 

.337 
Oil O 
.663.011 

it 
.989 
111 



Ordering the elements in the following way 

{(0, 0, 0), (0, 1, 0), (1, 0, 0), (1, 1, 0), (0, 0, 1), (0, 1, 1), (1, 0, 1), (1, 1, 1)} 



.451 .549 



.378 

.995 




Where, 






















.544 .456 
.337 



.622 


.998 






.663 





.005 


















.002 

.544 .456 

.337 

.113 .448 .439 

.011 












.663 


.989 



.113 .448 .439 
.011 .989 



. Observe that, 



.005 -.005 

.001 -.001 

.002 -.003 .005 

.002 -.002 



As a consequence, we obtain max\{Ti)ij — (T^)^! < .005, then the networks are 
.006-homeomorphiscs. The steady state of Ti is tti = (0, .01926, 0, .98074), and 
the steady state of is tt^ = (0, .01632,0, .98368). We can see that \tti — tt^I = 
maXi\Tri{i)-TT^{i)\ < .004. Meanwhile for T2 we have, tts = (0, 0, 0, .01632, 0, 0, 0, .98368). 
We can observe that with the new order in the elements of Z2^, the transition ma- 

Tii T12 



trix is an invariant submatrix of T2, that is T2 = 
example, Xi and <p{Xi) are e-similar. 



Ts 



In the above 
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To introduce the following theorem we use the transition matrices in this exam- 
ple, and we can observe the following calculation 



therefore max \{Ti) 



Tt - T/ = 



-.001467 


-.00232 




<.003. 

.000394 


.000161 




-.00136 
.00199 
-.00019 
.002639 



-.00044 
.002525 
.00156 
.002843 



.00006 


.00295 




.00011 


.00213 




.00277 
-.00199 
-.00243 
-.00263 



.00073 
-.00253 
-.00353 
-.002843 



and max\{T^),j - iT^)^j\ < .004. 

Theorem 2.2. If (p : Xi Xi is an e-homomorphism, then the transition matrices 
Ti and satisfy the condition 

maa;|(ri")„-(V)y|<e 

for all possible i and j, and all n > \. If the homomorphism is injective and 
e < 1, the steady state of Ti and the steady state of are close, that is satisfy 
ki - TT^I = maXt\TTi{i) - ■K^{i) 

Uo , therefore we have the fol- 



when /"(u^ 



Proof Since (Ti")y 
lowing: 

\p{uj\u))-p{cj,{u),cl,{f{u))\ = 
\p{u,f{u))p{f{u)j\u)) < 
fiu))\\piu, fiu)) - pi^u), Hfinm + 

Using this property, mathematical induction over n, and the cases giving for the 
possible values of the probabilities p{u, v), we can conclude that our aim holds. □ 

An e-homomorphism between two PRN determines a correspondence between 
the Markov Chains of these two networks. The following new definition is connected 
in some way with the idea of bisimulation for two TDMC (Time Discrete Markov 
Chains) studied and introduced in Two TDMC are equivalent by a bisimulation 
if they simulate the same problem with two different set of probabilities. Here, we 
introduce the concept of two similar TDMC 

Definition 2.3. The two TDMC of the same size n x n: {Ti, T^, Tf, . . .}, and 

{T2, Tf , Tf 7 • • ■} '^'"s e-similar or e-isomorphic if 

(^■ij)nxn satisfies 



(1) there exists an e G M enough small, such that T"' 
that \tij I < e, and X^iLi ^ij ~ Z*"" ™' 

(2) x{Ti^)ij = x{'^2^)i3! for all fn, where x is the characteristic function. 
That is, simulating the networks by these two TDMC are e-similar. 

In the above example, the TDMC generated by T and are .005-similar, and 
the networks simulated by them are .005-similar. 
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3. Construction of Probabilistic Regulatory Networks, examples 
Sum of two PRN 

Let Xi = {Xi,F= {fi)i^i,C) and X2 = (^2, G = igj)"U,D) be two PRN. The 
sum Xi®X2 = (X1UX2, F V G, C V L») is a PRN where 

(1) Xi(jX2 is the disjoint union of Xi and X2. 

(2) the function hij = {fi V gj) is defined by hij(x) — fi{x) x & Xi and 
hij (x) = gj (x) if X G X2 . 

(3) the probabihty ) = CiVdj, that is p{hij) = Ci if hij = fi or p{hij) = dj 
if h,j ^gj. 

If Ti and T2 are the transition matrices of Xi and X2 respectively. Then T = 
T \ 

„ ^ ) is the transition matrix of Xi (B X2. 

U J-2 J 

Example 3.1. 

An example of sum is the PRN obtained by summing the same PRN twice, 
X ® X. To make the disjoint union, we subindicate X with for the first X and 
with 1 for the second X. That is, the new set is 

XoUXi = {(0, 0, 0), (0, 1, 0), (0, 1, 0), (1, 1, 0)} 
U{(0, 0,1), (0,1,1), (0,1,1), (1,1,1)}. 

The digraph is: 

(1.1.0) (1,0,0) (0,0,0) 

(0,1,0) 

(1.1.1) (1,0,1) (0,0,1) 

(0,1,1) 

This is a way to construct a PRN over {0, 1}" using either one or two PRN over 
{0, l}"-\ since 2"-i + 2""! = 2". 

Superposition 

It is clear that a PRN is the superposition of several Finite dynamical Systems 
(FDS)^ni over the same set X with probabilities assigned to each FDS. Since each 
functions defined over a finite field can be wrote as a polynomial function, we will 
use this notation for functions over a finite field, li X = {0, 1} = Z2, the finite 
field of two elements, all the FDSs over X have one of the following state space, 
where fi{x) = x; /2(x) = 1; /3(x) = 0; /4(x) = x + 1, Vx £ X: 



Li 


L2 


L3 


u 
















i 


T 


n 


10 


1 


1 


1 
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If Pi denotes de probability assigned to Li, and Ti denotes its transition matrix, 
then the set of all PRN is described as follows. 

We denote by L1L2 the superposition of Li and L2, and similarly L1L3 is the 
superposition of Li and L3. The following state spaces are the superposition of two 
FDS with two elements: 





LIL3 


LIL4 


L2L3 


L2L4 


L3L4 


0^"^ 







OP' 





0^=* 


IP' ' 


-|-P3 ' 


P4. , 

Ti 




|P4-|-P2+P4 ' 


P4 

1 |P4+P3 


1 0^ 


1 0^^ 


1 


1 OP' 


1 0^=" 


1 



For example, with transition matrices 

r,. = T.+r. = (« -)r., = T. + T3=(;^ ) 

Product of two PRN 

Let Xi = {Xi,F= {fi)f^i,C) and X2 = (^2, G = {gj)jLj^,D) be two PRN. The 
product Xi X X2 ^ {Xi X X2, F X G,C A D) is a, PRN where 

(1) Xi X X2 is the cartesian product of Xi and X2. 

(2) the function hij = {fi,gj) is defined by 

hij {xi,X2) = {fi{xi),gj{x2)) 

for xi e Xi, and 2:2 € X2. 

(3) the probability is a function of Cj and dj, for example p(/i,j) = '''^^^ ■ 

Example 3.2. 

The product L1L2 x ^1X3 is the PRN with four states {(0,0), (0, 1), (1,0), (1, 1)} 
and four functions 

fii{x,y) = {x,y),fi3{x,y) = {x,0), 

f2i{x,y) = (l,y), f23{x,y) = (1,0). 
The state space is the following: 

L\L2 X L1L3 

P11+P13 Q (0,0) ^P'' (0, 1) OP'' 

P23+P21 I / IP21 
-I- 1/ P23 -I- 

(1,0) < — (1,1) QP'-'+P'' 

P13+P23 

The transition matrix is the following 

^ Pn + Pi3 P23 +P21 \ 

rp ^ P13 Pn P23 P21 
1 

\ JJ13 + P23 Pn + P21 J 
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3.1. Linear Probabilistic Regulatory Networks. A linear PRN is a superpo- 
sition of linear FDS. A linear FDS is a pair {X, f) where / is a linear function, and 
X is a vector space over a finite field. So, a linear PRN is a triple {X, {fi)^i,C), 
where X is a finite vector space, the functions fi-.X^X are linear functions, and 
C = {cj = The set X has cardinality a power of a prime number and each 

linear function is determined by its characteristic polynomial and the companion 
matrix. 

li X ~ — {0, 1, 2} is the field of integer modulo 3, then the linear functions 
are: fi{x) = x, f2{x) = 2x, and fsix) = for all x £ Z3. So, the linear PRN are 
the following: 

{/l)/2} {fljfs} {/2,/3} {/l,/2,/3} 

0l ^1 Ql Ql 



5P1 1 ^P2 2 OPi 1 Op, 2 1 2 Op, 1 ^p, 2 Op, 

If X = Z2 X Z2 is the vector space with 4 elements over the field Z2, then there 
are 4 linear FDS not isomorphics. In fact, using matrix, the possible characteristics 
polynomials p/(A) are: A^, + A + 1, A^ + A + 1. The companion matrices of 
these linear functions are: 



J ^ \ 1 J V 1 ; * V 1 1 

Then the FDS associated to this matrices are: 

Ai A2 
0(0,0) ^ (1,0) 0(0,0) ^ (1,0) 

T \ 

(0,1) (1,1) 0(0,1) ^ (1,1) 

^3 Ai 

0(0,0) 0(1,0) 0(0,0) (1,0) 

/ T 

0(0,1) 0(1,1) (0,1) ^ (1,1) 
The linear PRN with two functions are the following: 

Al,A2 Ai,A3 ^l'^4 

^(0,0)^(1,0) 6(0, 0)^(1, 0)8 ^MO, 0)1^(1,0) 

0(0, i)&(i, 1) o (0, 1) (1, 1) ^Y^^^^ 

A-2,A2. A2,Ai ^3,^4 

O^ (0, 0) (1, 0) op^* O^ (0, 0) ^P2 (1, 0) O^ (0, 0) (1, 0) 0P3 

Pi ^|P4 Pi y[P* 

O^ (0, 1) ^P^ (1, 1) OP' OP^ (0, 1) (1, 1) P3 (0, 1) ^P* (1, 1) 0P3 

Example 3.3. Looking the state spaces of the Ai,Aj for i ^ j PRN, we can 

conclude that two of them are not isomorphic. Working with the transition matrices, 
we can calculate when there exists an homomorphism between two of them. For 
example: the set Hom{AiA2; A1A3) = 

Example 3.4. 
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If Xi and X2 are two PRNs then the product Xi x X2 can be projected over 
each component. In the example of product of two PRN, we use the PRN L1L2 
and L1L3 to construct -L1-L2 x -L1-L3. The functions 4>i{x,y) = x and 4>2{x,y) = y, 
the usual projections, are e-homomorphism. In fact 

: Z2 X Z2 ^ Z2, and 0i o /n = /i o (/)i,and 0i o /21 = /2 o i/ii. 

(2) the second condition is satisfies too, since several arrows compact to one 
arrow in Z2, for example (/i'i(0, 1) = (^i(0,0) = 0, then the arrows compact in the 
following way 

{P11+P13Q (0,0) (0, 1) 0"''} 

Similarly 

{O^ (1,0) ^P^^+P^^ {1,1) OP^'+P^'} 101 
(0,0) (0.1) ] 

(1,0) (1,1) J 1 

(3) for the e-condition, the probabilities need to satisfy: 

bii-Pi|, \Pi3-pi\, \{pn + Pis) - Pi\ < e 

|(pi3 +P23) - 1|, |(Pii +P21) - 1| < e 
\P21 -P2I, \P23 -P2\, \{P21 +P23) -P2\ <e 

since pi + P2 = 1, and pi + = 1, for pij = piPj for example, we have \{pii + 

P13) = \{P21 +P23) -P2I = 0, \pii -pi\ = |P23 -P2I = blP2|, |Pl3 -Pl\ = \pI\, 

\P21 -P2\ = \pI\, |(pi3 +P23) - 1| = \Pi\, and |(pii + P21) - 1| = \P2\ = IpsI- So, 
max(pi,p2,P3) = e- 

We have the same conditions for 02. It is clear that always the projections are 
homomorphisms, but the applications of the e-condition depend on the particular 
case that we are working. It is clear that there are several ways to define inclusions 
i : L1L3 — > X LiL2- 

The following digraph shows the projections in general. 

Xi X X2 

^ TTl 7r2 

X\ X2 

The functions tt, : Xi x X2 — > Xi, are defined are follows: 7ri{xi,X2) = Xi, i = 1, 2. 
It is easy to check that the others two conditions are satisfied as in the above 
example. 

4. Invariant Subnetworks and Projections 

A subnetwork Y C X oi X — {X, F, C) is an invariant subnetwork or a sub- 
PRN of X if fi{u) G Y for all u € Y, and fi G F. Sub-PRNs are sections of a PRN, 
where there aren't arrows going out. The complete network X, and any cyclic state 
with probability 1, are sub-PRNs. An invariant subnetwork is irreducible if doesn't 
have a proper invariant subnetwork. An endomorphism is a projection if tt^ = n. 

Theorem 4.1. // there exists a projection from X to a subnetwork y then y is an 
invariant subnetwork of X. 
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Proof. Suppose that there exists a projection tt : X —>■ Y. If y S by definition 
of projection TT{y) = y, and /i(7r(y)) = ■JT{gj{y)). Therefore all arrows in the 
subnetwork Y are going inside Y, and the network is invariant. □ 

Example 4.2. 

The PRN X has two invariant subnetworks with projections 

■^i{x,y) = {x,0) and Tr2{x,y) = il,y). 



.67 

S ^0 



•67 , 

O (0,0) „ ^1 

1 .33 I ^ .33 I "^l ■* 
1 6(1,0) 



■67 , 

O (0,0) 



.46 



(0,1) o 



39 .68 

_ (1,0)^(1,1) 
X 

i 

S2 0^ (1,0)^ (1,1) o 



.68 



Checking the probabilities for m and 7r2, we have ei = .68; and 62 = .67. We can 
observe that X = Si x 82- 

Example 4.3. 

The subnetwork Xi = {{{x,y,l)}, F,C) is an invariant subnetwork of X = 

{{o,iy, F,c). 



000 



X 



.549 

.451 .995 
\ i 

010 

.378 

o 



.005 



.622 



001 

110 

Q.998 



■ 1 13, .4 56 

.544 

\ 
.002 



.448 

101 o 

.439 

I 

.337 

Oil 

.663.011 

it 
.989 

111 o 



00 



.1 13, .45 6 
.544 



Xi^Xi 



10 

.439 

i 

01 

.011.663 

Ti 
11 



.448 

o 



.337 

o 



Ordering the elements in the following way {(0, 0, 0), (0, 1,0), 
(1, 0, 0), (1, 1, 0), (0, 0, 1), (0, 1, 1), (1, 0, 1), (1, 1, 1)}, the matrix 
.544 .456 
.337 .663 
.113 .448 .439 
.011 .989 
Til T12 



is an invariant part of the transition matrix 
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Using the projection t: : X ^ Xi, 7T{x,y,z) — (x,y,l); and the isomorphism 
p{x,y, 1) — {x,y), the network X is projected over the network Xi. Checking the 
arrows the projection 7f is a .5-homomorphism. 

Example 4.4. The following PBN appears in flT, and has three suh-PBNs. 



1/ 
(100) 

T 



(010) 

(P2+P4) 



1 



P2 



(Oil) 

I 

P4 

i 

1 O (000) 



(001) 

P3 

/ 

Pi 



\1 

(110) 

T 

P1+P3 P2+P4 
i I 

(101) 

I 

P1+P3 

i 

(111) 



Xi = {(000)}, X2 = {(111)}, and X3 = {(100), (010), (110), (101), (111)} are sub 
networks. With adequate order the transition matrices Tx and are: 



T. 



X3 




1 

.5 




.5 








.5 











.5 
1 



and TTX; = (0,0,0,0,1) 



meanwhile Tx — 



Til T12 
Tx, 



, and Tlx = (.25, 0, 0, 0, 0, 0, 0, .75). 



5. The category of Probabilistic Regulatory Networks, and 
mathematical background 

Theorem 5.1. If (pi : Xi X2 is an ei-homomorphism, and (j)2 '■ X2 ^ X^ is 
another e2-homomorphism. Then (j) — 4>2 ° (pi '■ '^1 ~^ ^'i is an e-homomorphism. 
Therefore the Probabilistic Regulatory Networks with the homomorphisms of PRN 
form the category PRN. 

Proof. The Probabilistic Regulatory Networks with the PRN homomorphisms is a 
category if: the composition is an homomorphism, and satisfy the associativity law; 
and there exists an identity homomorphism for each PRN. 

(1) Let (pi : Xi ^ X2 be an ei-homomorphism, and let (p2 '. X2 X3 be an 
e2-homomorphism. If qt, gk and fj are functions in each PRN, and such that 
(pi o fj = gk o (pi and (p2 ° 9k — It ° <p2, then we will prove that: (p o fj = qt o <P- In 
fact, 

{(p2 O (pl) O fj =(j)20 {(pi O fj) = 

4>2 o (gk o (pl) = i(p2 ° gk) °(pi ^ 
(qt o (P2) °(pi = qt° (02 o (pl)- 

(2) We want to prove that 

X{tk{.(p{u),(t>{v))) > x(q(m,w)). 



PROBABILISTIC GENE REGULATORY NETWORKS, ISOMORPHISMS OF MARKOV CHAINS3 



Suppose that x(ci {u,v)) — 1. Then, since 0i is an homomorphism of PRN, we have 
that 

X{dj{(t>i{u),(l)i{v))) > x{ci{u,v)) 
which is 1. Since (p2 is an homomorphism of PRN, we obtain that 

>xicj{Mu),{Mv)) = 'i- 

Therefore we obtain that 

xitk{MMu)),MMvm = ^- 

Then the composition of two PRN-homomorphisms is an homomorphism. 

(3) To verify the third condition for e-homomorphism, we do the following. If 
p{4>{ui), (j){u2)) > 1, with ui, U2 G -'^i, then we need to prove that there exists an 
e such that 

\p{ui,U2) -p{(p{ui),(j){u2))\ < e. 

In fact: 

\p{ui,U2) - p{(t){ui),(l){U2))\ = \P{U1,U2) - p{(t)l{ui),<f>l{U2)) + 
p{(t)l{ui),(f)i{U2)) - p{(f>2{(t)l{ui)),(f)2{(t)l{U2)))\ 
< \p{ui,U2) -p(0l(ui),0l(U2))| + 
\p{(f)l{u-i_),(l)i{U2)) - p{4>2{Mui)),<t>2{'Pl{U2)))\ < ei + £2 

\p{ui,U2) - p{(j){ui), (/'(w2))| < ei + e2 

because and ^2 are e-homomorphisms. 

The associativity and identity laws are easily checked, therefore our claim holds, 
and PRN is a category. □ 

It is clear that, the PRN with the homomorphism between them form a category 
that we will denote VTZN . The category PRN is a subcategory of VTZJV, since 
an homomorphism is not always an homomorphism for some e G R enough small. 
But, if we don't include the condition for e to be enough small, the two categories 
are the same, because always an homomorphism is an e-homomorphism for some 
e G R. 

Theorem 5.2. XeiAfixAs = {X-ixX2,H,E) be a product of PRN Xi = {X-i,F,C) 

and X2 = {X2,G,D). If Si : X ^ Xi are two PRN-homomorphisms, then there 
exists an homomorphism d : X ^ Xi x X2, such that (piO S = 6i for i = 1,2. That 
is, the following diagram commutes 

Xi XX2 

<t>l 5 (t>2 

/ T \ 

X^^X^X2 

This homomorphism is unique. 
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Proof. The function 6 : X ^ Xi x X2 is defined as follows 5{x) = {Si{x), d2{x)), 
X G X. (5 is an homomorphism, in fact: 

(1) Let X = {X,L,P) be a PRN. Since 5i and S2 are homomorphism, for all 
function It & L there exist two functions fi £ F and gj G G, such that Siol^ = fio5\, 
and 62 °lt = gj ° ^2- Then for the function It there exists the function {fi,gj) that 
satisfies S o It = {fi, gj) o 5. 

{6olt){x) = 5{lt{x)) = {Si{lt{x)),52{lt{x))) = 

{MS,{x)),gj{S2{x))) = {{fi,gj)oS){x) 

(2) In order to prove x(eij((5(a;), (5(a;'))) > xiPhi^, x')), suppose x{Ph{x, x')) = 1. 
Then lt{x) = x', and 5{x') = 6{lt{x)) = {fi, gj){S{x)) by part (1). Therefore 
x{eij{S{x) , S{x'))) = 1, and our claim holds. 

It is easy to check that 6 = Si, in fact 

0i((5(a;)) = (t)i{Si{x),52{x)) = Si{x), 
for all xgX. □ 

If ^i, i = 1, 2, are ej-homomorphism then 

max\p{x,x') - p{(l)i{S{x)), ^i{S{x')))\ < ei. 

But 

\pix,x') -piSix),Six'))+piSix),S{x'))-p{MSix)),MSix'm < 
\p{x,x')-p{6{x).S{x'))\ + 
\p{6ix),S{x'))-p{Mm),MKx'm < ei. 

Therefore 

\p{x,x') -pi6ix),Six'))\ < ei - \pi6ix),6{x'))-p{Mm),MS{x'm 

\pix,x') - p{S{x),5{x'))\ < -Ii. 

Therefore 5 is an 6-homomorphism. So, the theorem holds for e-homomorphism. 

It is an inmicdiato consequence the following result, also is true for e-homomorpliisms. 

Theorem 5.3. Let Xi®X2 = {X1XX2, H,E) be a product of PRN Xi = {Xi,F,C) 
and X2 = {X2 ,G,D). If ji : Xi ^ X are two PRN-homomorphisms, then there 
exists an homomorphism 7 : Xi © X2 — » X , such that 7 o = 7^ for i = 1,2. That 
is, the following diagram commutes 

Xi®X2 

1-1 7 1'2 
/ i \ 

Xx^X^X2 

This homomorphism is unique. 

Theorem 5.4. All reducible PRN is either a product of its non trivial sub-PRN or 
a subnetwork of this product. 

Proof. It is trivial by definition of Product and sub-PRN. □ 
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6. Conclusions 

The intersection, and the union of two sub-PRN is a sub-PRN, therefore the class 
of sub-PRN of a particular PRN is a lattice. Reduction mappings described in [TT] 
and defined for PBN using the influence of a gene, for example Xn, on the predictor 
function /j'' to determine the selected predictor, can be extended to PRN. In order 
to extend this procedure to more than boolean functions, we use the polynomial 
description of genetic functions given in 0, the partial derivative is the usual in 
calculus and all the concepts in JI] can be using for PRN. Similarly our definition 
of projection, the reduction mappings are e-homomorphisms, and we can use for 
genes with more than two quantization, since this extension is not a trivial work 
we develop the theory and methods in 1 . 
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